Molecular vibrations-induced quantum beats in two-dimensional electronic 
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Quantum beats in nonlinear spectroscopy of molecular aggregates are often attributed to electronic phenomena 
of excitonic systems, while nuclear degrees of freedom are commonly included into models as overdamped 
oscillations of bath constituents responsible for dephasing. However, molecular systems are coupled to various 
high-frequency molecular vibrations, which can cause the spectral beats hardly distinguishable from those 
created by purely electronic coherences. Models containing damped, undamped and overdamped vibrational 
modes coupled to an electronic molecular transition are discussed in this paper in context of linear absorption 
and two-dimensional electronic spectroscopy. Analysis of different types of bath models demonstrates how 
do vibrations map onto two-dimensional spectra and how the damping strength of the coherent vibrational 
modes can be resolved from spectroscopic signals. 



I. INTRODUCTION 

Various ultrafast spectroscopy techniques are capable 
of probing molecular dynamics on timescales of nuclear 
motion and spectroscopic signals usually reveal effects 
caused by the coherent exciton dynamics, energy trans- 
fer, vibrational motions, etc. perceptible even for simple 
systems^"— . For example, coherent quantum beats ob- 
tained at broadband excitation conditions are displayed 
in spectra of small molecular complexes^ —, thus, reflect- 
ing an interplay between the absorbed energy transfer 
and relaxation. 

The two-dimensional (2D) electronic spectroscopy 
(ES) displays highly resolved decoherence dynam- 
ics in molecular aggregates induced by the system- 
bath interaction. Problem of assignment of the de- 
tected slowly-decaying spectral beats to either electronic 
coherences 3 -^ - — or molecular vibration s 14 ! 15 is under in- 
tense discussions up to now. However, it is still disputed 
whether such detachment is feasible at all. 

Molecules and their environment are characterized 
by a wide variety of high-frequency vibrational modes, 
some of them strongly-coupled to the electronic ex- 
citations. Unambiguous experimental evidences that 
the high-frequency vibrational wavepacket motion causes 
the spectral beats in 2D spectra were reported only 
recently^—. However, it is not yet clearly established 
how these vibrational resonances map onto 2D spectra 
and how the excited vibrational wavepacket evolution af- 
fects the spectral dynamics. 

Systems where the high-frequency vibrations are mixed 
with electronic transitions are currently of a special in- 
terest. Great efforts are made to effectively and univer- 
sally identify the underlying mechanisms^ —. Analysis 
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of rephasing and nonrephasing signals of the 2D ES sep- 
arately could help solving this issue in some cases, but 
usually the strong spectral congestion leaves ambiguity 
of the result due to the overlap of many oscillating fea- 
tures. In case when the limited-bandwidth laser pulses 
cannot cover the whole absorption spectrum even if the 
vibrational features are well-resolved, the intramolecu- 
lar vibrations can still be evaluated by analyzing the 
phase relationships of beating patterns in the rephas- 
ing and nonrephasing signals, peak ellipticity or rota- 
tion of the nodal line of the imaginary part of the total 
2D spectrum ^ 19 ! 21 . However, vibrational coherences can 
be created in the electronic ground state manifold re- 
sulting in beats of both diagonal and off-diagonal peaks, 
thus making the whole time-resolved dynamics even more 
entangled— »21. 

Non-decaying vibrations-related quantum beats on top 
of the vertical electronic transitions are obtained in sim- 
ulated time-resolved spectra when molecular vibrations 
are included as coherent harmonic modes of the bath 
coupled to the system 4 ^ 22 ! 23 . This problem sometimes is 
formulated in pure quantum mechanical terminology by 
using so-called polaron transformation 24,25 . Such treat- 
ment is exact when non-interacting molecules are con- 
sidered; however, in molecular aggregates vibrations are 
coupled to the electronic degrees of freedom in a non- 
trivial way. More sophisticated theoretical methods, 
for example, hierarchical equations of motion or time- 
dependent density matrix renormalization group tech- 
nique, might resolve the dynamics of such mixed vibronic 
systems 2 ^"—. As the vibrational coherence decay rate is 
usually smaller compared to those of the electronic co- 
herences, this might also trigger the reconsideration of 
the nature of long-lived coherences observed in many sys- 
tems. 

In this paper, we theoretically analyze the influence of 
high-frequency damped vibrational motion of molecules 
on the spectroscopic observables and make comparison 
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with the undamped description. In the microscopic 
model such treatment translates into a more sophisti- 
cated expression of bath spectral density function. Gen- 
eral bath damping mechanisms are also the problem un- 
der consideration; the model of a quantum overdamped 
bath is introduced and its influence upon the lineshape 
formation in the linear absorption and 2D spectrum is 
discussed. Strongly-overlapping spectra of simple elec- 
tronic and vibrational systems are examined in the con- 
text of various bath models. 



II. GENERIC MOLECULAR COMPLEXES 

Electric fields interacting with molecular systems can 
affect various system degrees of freedom ranging from 
electronic high-energy states to low-frequency molecular 
vibrations. It is natural to simplify the description by 
taking into consideration only the most relevant types of 
dynamics. One of the main simplification is achieved by 
neglecting the non-resonant transitions with respect to 
the optical field under consideration: in the case of the 
resonant excitation with the optical laser field a molecule 
is described as a system consisting of two electronic en- 
ergy levels, while a supramolecular complex - as an ag- 
gregate of coupled two- level electronic systems. Another 
standard assumption when describing electronic quan- 
tum levels of a molecule is the separation of electronic 
and nuclear degrees of freedom (the Born-Oppenheimer 
approximation). It is assumed that the electronic wave- 
function rigorously follows the nuclear configuration and 
the energies of electronic levels become parametric func- 
tions of nuclear degrees of freedom. The equilibrium nu- 
clear configuration which is usually different for each elec- 
tronic quantum state is defined as the energy minimum 
with respect to the nucle i 29 ! 30 . 



A. Absorption of an isolated molecule 

Using the assumptions made above, electronic exci- 
tations in molecules are approximated as the Franck- 
Condon transitions. It is considered that the (slow) nu- 
clear degrees of freedom remain frozen during the elec- 
tronic transition and the fixed molecular nuclear configu- 
ration adjusted to the ground electronic state emerges in 
non-equilibrium conditions with respect to the new elec- 
tronic state. As the electronic potential energy surface 
is assumed to be parabolic depending on the nuclear dis- 
placement in the vicinity of the equilibrium, a displaced 
(harmonic) oscillator model can be used in the descrip- 
tion of the electronic transition (Fig. [I]). In one dimension 
the electronic potential of the ground state is V g (q) = 
huooq 2 /2 and the displaced electronic excited state is de- 
scribed by the potential V e (q) = u eg + fowo(q — d) 2 /2. 
Here ujq is the vibrational frequency, uj eg is the energy 
gap between the minima of two potentials and d is a 
dimensionless displacement parameter determining the 




Fig. 1. Model system of displaced harmonic oscillators and 
Franck- Condon overlap factors for the corresponding m — n 
transitions. 

strength of the electron coupling with molecular vibra- 
tions. The Huang- Rhys factor s = \d 2 is widely used to 
qualify the coupling strength. 

Quantum mechanical problem of the isolated harmonic 
potentials has exact solutions. They constitute an infi- 
nite set of wavefunctions ip m with quantum number m — 
. . . oo and corresponding energies E m = ftwo(m + 1/2) 
with respect to the bottom of the corresponding potential 
surface. Transitions between the sub-states of the elec- 
tronic ground state and the ones of the electronic excited 
state determine the vibronic progression in the absorp- 
tion spectrum (Fig. [T]). The intensity of each vibronic 
peak is scaled by the overlap of vibrational wavefunc- 
tions - the Franck- Condon overlap - in the ground and 
excited state potentials. Following these assumptions the 
absorption spectrum of a single isolated Franck- Condon 
molecule can be given by 

CO 

< b C »OCW e~^\Fmn\ 2 (1) 

m,n=0 

POO 

xRe / ^j- e ^-uj eg )t-\ujQ{n-rn)t-^t 
JO 

where the line-broadening parameter 7 is introduced phe- 
nomenologically. F mn is the Franck- Condon wavefunc- 
tion overlap integral for the m to n transition 

Fmn = (m\ exp ^--^= (a 1 " - a)^j \n) (2) 

given in terms of bosonic (vibrational) creation a) and 
annihilation a operator ^ 29 ! 30 . 

B. Molecule coupled to the phonon bath 

In order to capture both the vibrational-type set of 
states of the energy spectrum and the spectral broad- 
ening in an unified model a more general description of 
a molecule is needed. The model discussed above con- 
sidered the composite electronic+vibrational system as a 
closed system. Consequentially, the line-broadening pa- 
rameter 7 was included phenomenologically and without 
more detailed physical insight. 
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The total Hamiltonian of a single molecule coupled to 
the set of phonon modes (an open quantum system) is 
given by 

A = 0\g){g\ + (u eg + \)\e)(e\ (3) 
-^hjjj^ ] y/s^ {a ] a + a a ) \e){e\ 



Here the ground state \g) is the reference zero-energy 
state, the energy of the excited state \e) is additionally 

shifted by the reorganization energy A = J2 a s & •> 
while the rest of vibrational degrees of freedom are repre- 
sented by a set {a} of harmonic oscillators. The system- 
bath coupling represented by the third term in Eq. (j3]) 
introduces fluctuations into the energy of the electronic 
excited state due to low-frequency bath modes at fixed 
temperature. The most convenient form to describe such 
fluctuations is the time correlation function: 



C(t) =Tr B {Oe(t)Oe(0)peq} 



Here 



exp 




(4) 



(5) 



is the normalized thermally equilibrated canonical den- 
sity operator of the bath (/3 _1 = fceT), and 



Q e (t) = OaW + Oa(*)) 



(6) 



is the fluctuating collective bath coordinate in the Heisen- 
berg representation. The exact form of the correla- 
tion function describes the high-frequency modes as well. 
Trace in Eq. (|4]) amounts to averaging over the infinite 
set {a} of harmonic oscillators. Inclusion of eqs. ([5j) and 
(|6|) into the trace gives 



C(t) = ^h 2 {oj^y s a (coth 

a \ 



(«) 



C'(oj) and C"(oj) are related by the fluctuation-dissipation 
theorem 



C'(w) = coth(phu]/2)C"(w). 



(10) 



To describe the optical properties of the molecule the 
system polarization operator P = \g) (e\ + \e) (g\ is as- 
sumed (the transition dipole strength is taken as unity). 
It can be shown that the quantum correlation functions of 
the polarization operator representing the spectroscopic 
observables are exactly given in terms of the second- 
order cumulant expansion with respect to the vibrational 
modest. Therefore, it is exact for our molecule and we 
can calculate the shapes of electronic transition bands 
and the bath-induced time dependence of the spectrum 
by using the formalism of the lineshape functions. The 
latter directly comes from the perturbative second-order 
cumulant expansion of the system density operator prop- 
agation which allows to describe various types of vibra- 
tional baths and include these effects explicitly 31 . The 
absorption coefficient is then given by 



JO 



(ii) 



Here dimensionless lineshape function g(t) is an integral 
transformation of the correlation function C(t) of system- 
bath fluctuations, or for its Fourier transform (Eq. (|5J)), 



g(t) 



2^ / 



C{u>) 



[exp(-icjt) + loot - 1] . (12) 



Assuming that the system is coupled to a continuous 
spectrum of bath frequencies, the correlation function 
Eq. © can be calculated for a predefined distribution. 
We introduce the density of the Huang-Rhys parame- 
ter as a function of vibronic frequency, {s a } — >> s(u)duj, 
which may have a peak at some dominant normal-mode 
frequency ujq. Then, from eqs. (0) and ((9|) the spectral 
density is obtained as 



C"{u) = TrfrV [s(lj) - s(-u)} . 



(13) 



■ coscJq J t — lsmc^Q J t 



(7) 

which is the well-known form of the two-point correlation 
function of generalized bath coordinates 30 . The Fourier 
transform of the correlation function is a real function 



For instance, if we consider a vibrational damped mode it 
will be represented by a broad peak in the spectral den- 
sity of the system-bath coupling. In case of the Gaussian- 
type coupling 



1 jazz 
e 2 



CM 



27T7 



dte iu;t C(t) =C'(w)+C"(u), 



(8) 



the spectral density function is 



where C'(uo) and C"(uo) are even and odd functions of 
oj. C"(uo) is the temperature-independent function and is 
denoted as the spectral density: 



Cq(uj) = A/i- 



a/27 



2 7 2 



while the Lorentzian-type coupling 



/'OC 

C"(w) = -2 / sinojtlmC(t)dt. 
Jo 



(9) 



1 



1 



7T (UJ - UJq) 2 + 7 2 



(14) 



(15) 



(16) 



leads to 



(w 2 -w2_ 7 2) 2 + 4w 2 7 2 



(17) 



Here A and \i are just some scaling constants here in- 
serted for convenience: the expression of \i is later chosen 
in order for A to be equal to the reorganization energy 

By this point, no approximations were applied - equa- 
tions (fT5j) and (p~7j) are consistent with the fluctuation- 
dissipation theorem. The inverse Fourier transform of 
them would give the time correlation functions which 
in case of Gaussian and Lorentzian couplings decay as 
e -7 1 / 2 and e _7t , respectively. The spectral densities in- 
clude the damping parameter 7 and vibrational frequency 
cjo- By taking various limits with respect to these param- 
eters, different damping regimes can be achieved repre- 
senting different conditions of the bath and both spectral 
densities (eqs. JTSJ) and (|T7|)) can be used in numerical 
simulations. 

We further analyze different regimes of undamped , 
damped and overdamped vibrational motion considering 
the Lorentzian-type coupling since it is related to the 
exponentially-decaying time correlation function, which 
is a natural decay pattern in most of physical applica- 
tions. When a single bath mode is assumed, i. e. when 
the spectral density function is obtained as a Fourier 
transform of a single term of Eq. ([7j) , the spectral density 
is given by 

C^(u) = 7tsu 2 [5(u - uj ) - S(uj + ljq)] . (18) 

The reorganization energy in this case is sujo and the 
lineshape function is 



ffu(t) 



coth^a 



cos coot) + i (sinc^o^ — ^0^) 



(19) 

This vibrational mode leads to non-decaying, undamped, 
vibrational motion. Such spectral density given by the 
^-functions is not realistic because dissipation, which ac- 
companies any vibrational motion due to the molecular 
interaction with its environment, is neglected. Damp- 
ing induces the decay of the correlation function over 
time and the corresponding spectral density should have 
a finite-width smooth peak. Such damped regime is 
achieved by taking 7 < ujq in Eq. (fT7|) and leads to 



c» = 



2V2A, 



+ 2 7 2 w 2 



(20) 



This spectral density function with reorganization energy 
A has a peak at ujq and the peak width is defined by the 
damping strength 7, differently from Eq. (fT8|) . where the 
peak is the function (Fig. [2^). 

The opposite, overdamped , regime is usually repre- 
sented by the spectral density of a Brownian oscilla- 
tor, which represents the fluctuations of the Ornstein- 
Uhlenbeck process^ 2 -. This regime can be postulated 




(b) 





7 

0.05^0 






0.25^0 






0.5c^ 






0.75^0 












1.5u; 





UJ — UJ t 



eg 




Fig. 2. (a) Spectral densities of undamped (red dashed line) 
and damped (black solid line) baths, (b) Absorption spec- 
tra of the displaced harmonic oscillator (s = 0.3) for different 
damping strengths 7/cJo = 0.05 — 1.5 of the damped bath 
(black solid lines) and for the undamped bath (red dashed 
line). The widths of three main peaks are displayed as the 
functions of damping strengths in (c), where the third peak 
becomes unresolvable for 7 > 0.5cjo- The real and imagi- 
nary parts of the lineshape functions corresponding to the 
different vibrational bath models are plotted in (d) and (e), 
respectively. 



semi-classically, i. e. by using the exponentially decaying 
classical correlation function C c \(t) = 2A/cBTexp(— j\t\). 
Its Fourier transform C c \(uo) = A^\k^T(uo 2 + 7 2 ) -1 now 
represents the even (classical) part C'(uo) of the total 
quantum correlation function. Since the classical trajec- 
tory reflects the high temperature limit, we have C"(uo) = 
{yS) and obtain the spectral density representing 
the classical Brownian particle 



.30 



2A7CJ 
uj 2 + 7 2 



(21) 



with reorganization energy A and the relaxation rate 7. 
The full quantum correlation function in the frequency 
domain can now be constructed by the direct application 
of the fluctuation-dissipation relation Eq. (fTU|) . There- 
fore, we denote bath described by such spectral density 
as the overdamped semi-classical bath. At the high- 
temperature limit the corresponding lineshape function 
is obtained as 



=(*) 



2 



(e-T* + 7 i-l). (22) 



Evidently, the overdamped bath can also be obtained 
from Eq. (fT7|) . In the limit of 7 3> loq it yields rather 
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different expression 



4Ac^7 3 



(uj 2 +7 2 ) 2 ' 



(23) 



We denote this regime as the quantum overdamped . The 
reorganization energies are equal to A for both types 
of the overdamped spectral density. The semi-classical 
spectral density function is equal to A at its maximum 
at uj = 7, while the quantum function has its maxi- 



mum value C"-q 



(^peak) A w 1.3A, where cj pea k 



a/3/37 ~ O.587 ( see Fig. [3^). The quantum overdamped 
spectral density decays faster at infinity, therefore, it is 
more suitable for numerical applications. At the high- 
temperature limit the corresponding lineshape function 
is given by 



go 



-qW = Z ( £ - 1 ) ( e " 7 V + 2e"^ - 2 + 7*) (24) 



A /_2_ 

7 



e 7£ + 7 £- 



III. SPECTROSCOPIC OBSERVABLES 

The absorption signals of the two-level system coupled 
to the bath are calculated according to Eq. (fTT|) and us- 
ing lineshape functions corresponding to different damp- 
ing regimes (eqs. (fT5|) . (|2Q|) . (f2TT) or ([23])). The included 
overdamped modes lead to the Lorentzian or Gaussian 
lineshapes of spectral peaks for fast , 2Xf3~ 1 j~ 2 ^> 1, or 
slow , 2A/3 -1 7 -2 <C 1, decay regimes, respectively. Fast- 
decaying modes of molecular vibrations will result in ho- 
mogeneous broadening, while the strong coupling to the 
high-frequency vibrations will result in vibrational pro- 
gression in the absorption spectrum. 

In the case when the molecule is coupled to a single 
high-frequency mode and a continuum of low-frequency 
damping modes the spectral density consists of two parts: 



c» = c>) + c>M- 



(25) 



Here the (fast) C"(uo) mode corresponds to the 
semi-classical C"_ sc (lj) or quantum C"_ q (u) overdamped 
bath. The second term, C" ib (cj), represents spectral den- 
sity of (slow) molecular vibrations ( undamped C"{uj) or 
damped C^(uj)). The total lineshape function, which is 
a linear transformation of the spectral density, then con- 
tains two parts as well. 

The spectral density and lineshape function formalism 
presented above is extended to the photon echo signal cal- 
culations for the 2D ES 33 . The response function theory 
denotes the response function S(t 3 , t 2 , t\) as a sum of dif- 
ferent pathways of system density operator propagations. 
In this study we assume the impulsive limit (laser pulses 
of infinitesimally short duration), where the 2D spectrum 
is defined by the two-dimensional Fourier transform of 



the response function over variables t\ and t 3 \ these then 
have the meanings of the time delays between the first 
and the second laser pulses, and the third laser pulse 
and the detection time, respectively 3 ^. t 2 is the fixed 
time delay between the second and the third interaction 
and is denoted as the population time. The rephasing k\ 
(Fig. [4]) and nonrephasing k\\ (Fig. [5]) signals are distin- 
guished as the ordering of the interaction with the first 
two pulses is changed. The total signal is denoted as the 
sum of k\ and k\\ parts. Evolution of its diagonal and 
off-diagonal peaks closely follows the dynamics of corre- 
sponding populations and quantum coherences of the sys- 
tem density operator. In 2D spectra, fast-decaying modes 
of molecular vibrations will result in anti-diagonal peak 
(homogeneous) broadening, while the strong coupling to 
the high-frequency vibrations will result in multiple os- 
cillating diagonal and off-diagonal peaks reflecting the 
configuration of vibrational state o 35 i 36 . 

For a single two-level system the complete expressions 
for the rephasing and nonrephasing contributions of 2D 
spectrum can be given in terms of the lineshape func- 
tions: 



Ski(t 3 ,t 2 ,ti) 



1 



-iu>o(*3 — *l) 



(26) 



xe 



xe 



7 m (t 1 +t 2 )-9*(t 1 )+g*(t 1 -\-t 2 +t 3 ) 
Re[g(t 2 +t 3 )+g(t 2 )-g(t 3 )] 



and 



x cos {Im [g(t 2 + t 3 ) + g(t 2 ) + g{t 3 )]} 



(27) 



xe 



xe 



g(ti+t 2 )-g(ti)-g(t 1 +t 2 +ta) 
Re[g(t 2 +t 3 )+g(t 2 )-g(t 3 )] 



x cos{Im [g(t 2 + t 3 )+g(t 2 ) - g(t 3 )]} , 

respectively. It should be noted that these expres- 
sions are exact for the harmonic bath and carry no 
approximation^. 

For a multi-chromophoric system the response function 
is a sum of many contributions, which can be classified to 
the excited state absorption, excited state emission and 
ground state bleaching pathways. Additionally, the path- 
ways responsible for population transfer are present^. 
Here in eqs. (f26|) and (|27|) . there are only two pathways; 
the excited state absorption contribution producing neg- 
ative peaks in the 2D spectrum is not available due to 
absence of the third electronic state. 



A. Quantum vs. semi-classical overdamped bath 

The differences in the two-level system absorption 
lineshapes obtained by using solely the quantum or 
semi-classical overdamped bath models (C"(oj) = C"(u)) 
are minor (Fig. [3Jd) despite of a substantial difference in 
the shapes of spectral densities (Fig. [3^). The position 
of the absorption spectra are mainly determined by the 
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Fig. 3. (a) Overdamped quantum (red dashed line) and semi- 
classical (black continuous line) spectral density functions of 
an two-level system and corresponding (b) absorption spectra, 
(c-d) real and imaginary parts of the lineshape functions. 



imaginary part of the corresponding lineshape function 
(Fig. [3]i): for the semi-classical bath, it is lower by A/7 
compared to the quantum bath at jt ^> 1. The slopes of 
the imaginary parts are, however, both equal. The slopes 
of the real parts of the lineshape functions, which deter- 
mine the absorption linewidth, are different (Fig. [3b) . 
For g Q - sc (t) the function slope at yt ^> 1 is while for 

g Q - q (t) is at least twice larger, l^ + ^f 2 . Such differences 
result in larger lineshape broadening in the case of the 
quantum bath. 

The shapes of functions describing the spectral peak 
dynamics of the rephasing 2D spectra obtained by 
using the quantum model of the overdamped bath 
(Eq. ([24]) ) coincide with the semi-classical bath simula- 
tions (Eq. ([22]) ) at short population times, as demon- 
strated in Fig. [4)3 and Fig. [Hi. However, in the non- 
rephasing (Fig. [5b and Fig. [5]i) spectrum the quantum 
overdamped bath changes both the amplitudes of oscil- 
lations and the averages of peak intensities. In the 2D 
spectra the quantum overdamped bath model produces 
lineshapes with slightly smaller broadening. By adjust- 
ing the parameters both models can be tuned to reflect 
experimental broadenings. The quantum model has a 
short tail in the spectral density, which is advantageous 
for numerical simulations. 



B. Undamped vs. damped harmonic bath 

The properties of undamped and damped harmonic 
baths are demonstrated by assuming the composite spec- 
tral density given by Eq. ([25]) and setting the Huang- 
Rhys factor of the coherent mode to s = 0.3. In absorp- 
tion, such system demonstrates three well-resolved peaks 
of vibrational progression at frequencies uo = cj eg , o) eg +ujo 
and o) eg + 2uo (Fig. [2b)- 

The case when the composite spectral density is 
C"_ sc (oj) + C^(cj) corresponds to the overdamped semi- 



classical bath with undamped vibrations. This is the 
commonly used approach and is discussed in detail in the 
literature ^ 19 ! 37 . The produced vibrational progression in 
the absorption spectrum is drawn as the red dashed line 
in Fig. [2}d. Unsurprisingly, the absorption signal calcu- 
lated from eqs. (pQ) and ([TT]) are then identical. 

The effects caused by damping of vibrations are often 
ignored, however they may raise significant implications 
in the nonlinear spectra. In the linear spectra, compari- 
son of the result obtained by assuming the spectral den- 
sity of Co-scM + C d M ( black solid lines in Fi £- Efc>) to 
the undamped vibrations reveals an almost perfect cor- 
respondence for 7 = 0.05cjo? which illustrates a smooth 
transition from one model to another. It is also notice- 
able that the peaks in the progression broaden gradually 
when the damping strength 7 is increased. However, the 
broadening is not uniform - peaks that are at higher 
energies are broadened more. This is evident by evalu- 
ating peaks width dependence in the damping strength 
(Fig. Eb), obtained by fitting the spectra with multiple 
Lorentzian functions. As a result, the total spectral line- 
shape becomes asymmetric: higher-energy shoulder of 
the vibrational progression is reduced due to additional 
broadening. Such pattern of broadening is consistent 
with the lifetime-induced decay of vibrational coherences 
at higher excitation energies. 

The peak lineshapes in the 2D spectra obtained by 
using the damped case show slightly larger broadenings 
compared to those of the undamped vibrations both for 
rephasing and nonrephasing signals (figs. [4] and [5]). The 
non-uniform broadening, as follows from the absorption 
simulations, shows up in the 2D spectrum, as well. For 
the vibrational bath we separately study two cases: the 
damped vibrations (7 = cjo/4, Fig. [4H-b and Fig. [5^- 
b) and the undamped vibrations (7—^0, Fig. [4]>d and 
Fig. [5t-d). In both cases the lineshape function repre- 
senting the semi-classical overdamped bath is added to 
the vibrational part. The reorganization energy in all 
cases is the same (2sco>o). 

The vibrational wavepacket (coherence) beats show up 
as temporal oscillations in the 2D spectrum. Decay of 
the coherences in case of the damped vibrations results 
in decay of cross-peaks (compare spectra for t 2 = 27rJq~ 1 
in Fig. [4^ and Fig. (4t): since all spectra are normal- 
ized to the maximum of the rephasing signal at £2 = 0, 
the intensities of the cross-peaks and the upper diago- 
nal peak are evidently lower than those of the main peak 
at (lj 1,0)3) = (uj eg ,uj eg ). The main differences in oscilla- 
tory dynamics of peaks in the 2D spectra are due to the 
damping-induced decay of coherences in the case of the 
damped vibrations. 

Considering the phases of oscillations corresponding to 
different peaks, it is clearly observed that the model of 
vibrational bath does not change the phase relationships 
of peak oscillations. In general, the phases of peak oscil- 
lations are not the same for all peaks in rephasing and 
nonrephasing signals. For example, in the rephasing sig- 
nal lower diagonal peak 'Dl' oscillates as the negative co- 
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Fig. 4. 2D rephasing spectra at £2 = 0, 27rj C) ~ 1 and Sttcj^" 1 
with damped (a) and undamped (c) vibrations on top of the 
semi-classical bath. Intensities of the peaks, indicated by 'Dl' 
(diagonal lower), 'Ch' (cross-peak higher), etc., are depicted 
as functions of population time £2 in case of damped vibra- 
tions (b) and undamped vibrations (d) with semi-classical 
(solid lines) and quantum (dashed lines) overdamped bath. 
All values are normalized to the maximum of the rephasing 
spectrum at £2 = 0. 



7T 6tt 8tt 
u t 



, 0.3 








; 0.2 

f 








0.1 






« r \ w a 


/o.c 


/erf 






J-°-J 










) 2-K 47T 


67r Sir iOn 



UJ t 



Fig. 5. 2D nonrephasing spectra at £2 = 0, 27rcj ( ^ 1 and SttUq 1 
with damped (a) and undamped (c) vibrations on top of the 
semi-classical bath; the intensities of corresponding peaks in 
case of damped vibrations (b) and undamped vibrations (d) 
with semi-classical (solid lines) and quantum (dashed lines) 
overdamped bath. All values are normalized to the maximum 
of the corresponding rephasing spectrum at £2 = (in Fig. [4j). 



IV. DISCUSSION 



sine function with the population time £2, while all other 
peaks oscillate with opposite phase (Fig.[4j3 and Fig. [Hi); 
in the nonrephasing signal (Fig. [5b and Fig. [5]i) both 
cross-peaks oscillate as the negative cosine, while both 
diagonal peaks - with the opposite phase. Such phase 
relationships are typical for oscillations of vibrational co- 
herences and are known as being strongly dependent on 
the Huang- Rhys factor; this property ought to be em- 
ployed in distinguishing vibrational and electronic coher- 
ences in spectral. 



We considered the smooth transition from undamped 
through damped to overdamped vibrations and their pos- 
sible manifestation in the 2D ES signals. For small values 
of damping parameters (7 < A) the undamped vibra- 
tional model can be used, but for 7 w A, effects of non- 
uniform peak broadening in absorption and the damping 
strength-dependent decay of coherences in the 2D sig- 
nal become significant. We also introduced the quantum 
overdamped spectral density model which has a conve- 
nient decay character with frequency and can be used 
to model broadening effects similar to the semi-classical 
model of overdamped Brownian motion. 
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Fig. 6. Absorption spectra of the monomer with damped 
(solid black line) and undamped (dashed red line) vibrations; 
absorption spectrum of a dimer (green dash dotted line). 



Understanding of the damping-related effects on the 
spectra is important when considering the realistic exper- 
imental measurements, where the mixture of electronic 
and vibrational quantum beats is expected to be ob- 
served. 

To demonstrate the relations we turn to simulations 
and comparisons of model systems that have strong spec- 
tral congestion and underlying system properties may be 
obscure. To make a constructive parallel to the typical 
experimental results, properties of total spectra are high- 
lighted rather than rephasing and nonrephasing spec- 
tra separately. We compare monomers coupled to (un- 
damped and damped) vibronic resonances of frequency 
ujq. The results are compared to the elementary molecu- 
lar aggregate - a molecular homodimer with resonant 
coupling | J | = ^ (Fig. [6]). The parameters of the 
monomer are the same as in Sec. IIIIB1 The 2D signals 
of the homodimer are simulated using the conventional 
response function theory for molecular excitons as de- 
scribed in Ref, 5 . The site energies of the dimer are taken 
to be e = uj eg + t;UJo and the inter-dipole angle (j) = |7r. 
Electronic transition to two excitonic states at uu eg and 
uj eg + ujo are then possible. The higher energy peak is 
more intensive than the lower one by factor ]~ cos t — 3. 

J 1 + COS0 

For such H-type dimer— the excited state absorption con- 
tribution will be much stronger in the lower cross-peak, 
making it negative in the 2D spectrum (Fig. [7^). The 
inter-dipole angle also determines the amplitudes of co- 
herence oscillations. Reducing the angle would increase 
the amplitude of coherence oscillations, but decrease the 
excited state absorption influence on the lower cross- 
peak. The spectral density corresponding to the over- 
damped Brownian oscillator (overdamped semi-classical 
bath, Eq. ([22]) ) is used for the spectral broadening. Its 
parameters are A = |cjo and 7 = ^ujo- 

The total spectra of the monomer with undamped vi- 
brations (Fig. [7)3) show a typical non-decaying 'breath- 
ing' behavior of the main pea k 16 i 17 . It can be charac- 



terized as the anticorrelation of diagonal and antidiago- 
nal widths of the main peak. Such change of the peak 
shape is created by the superposition of k\ and fen sig- 
nals, which have phase- inverted dependence of the peak 
intensity on £2, as can be seen in Fig. [7^ and Fig. [7f. The 
same relation is observed in simulations of systems with 
well-resolved peaks as shown in Fig. [4] and Fig. [5l 

As it was already mentioned, the phase relationships 
and oscillation amplitudes of the k\ and fen signals are 
Huang- Rhys factor-dependent, therefore, both correla- 
tion and anticorrelation of the diagonal peak- width and 
intensity is possible for different systems^. However, 
the anticorrelation of the diagonal width and intensity 
of any peak in the 2D spectrum is not available for the 
electronic dimer—. As follows from our calculations of a 
two-level system the amplitude is higher for the fei signal 
and, therefore, the correlation of the diagonal width and 
intensity of the main peak is observed both for damped 
and undamped vibrations. Thus, spectral dynamics of 
the vibrational and electronic-only systems are very sim- 
ilar by this aspect. 

The monomer with damped vibrations being analyzed 
here is an illustration of the transition from damped to 
overdamped vibrations. Due to damping, the main peak 
experiences only a few 'breaths' and the round shape is 
established at > 2.57tcj ( ^ 1 at least (Fig. [Tfc). Increas- 
ing the damping strength would destroy all vibrational 
coherences and it would not be possible to spectrally dis- 
tinguish it from a general two-level system coupled to 
an overdamped bath. However, for a molecular dimer 
one can expect it to demonstrate similar spectral char- 
acteristics compared to the monomer with damped vi- 
brations. Indeed, this is true for the 2D spectra pre- 
sented in Fig. [T^l. Due to the excited state absorption, 
a negative region below the main peak appears for short 
population times. The oscillation pattern of the peak 
maximum (green dot-dashed line in Fig.[7]i-f) closely re- 
sembles the dependence of the monomer with damped vi- 
brations (black solid line), as well. The decay rates of the 
signals are very similar here, but the established values 
of peak intensities at long population times are different. 
This can be explained by the additional dephasing mech- 
anisms (exciton transfer and lifetime decay) available for 
an excitonic system during £2. 



V. CONCLUSIONS 

Consideration of high-frequency damped molecular vi- 
brations instead of undamped is a more realistic descrip- 
tion since it includes dephasing and dissipation that vi- 
bronic motion experiences in solvent. This correction 
induces non-uniform peak broadening, changes the po- 
sition of peaks within the vibrational progression in the 
absorption and results in decay of coherences in the 2D 
spectrum. The significance of the effect depends on the 
width of the spectral density function around the reso- 
nance. Therefore, having the coherence decay rate and 
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Fig. 7. Real parts of total (ki + kn) 2D spectra of a dimer (a) and a two- level system with undamped (b) and damped (c) 
vibrations at population times £2 = 0...4.597rcj^" 1 . Plots of intensities of peaks, indicated by 'X' of total (d), k\ (e) and ku 
(f) signals of the monomer with damped vibrations (solid black lines), undamped vibrations (dashed red lines) and molecular 
dimer (dot dashed green lines) are depicted. Spectra are normalized to the values of intensity maxima at ti — for each system 
separately. 



the peak broadening quantitatively evaluated from the 
experiment, one would be able to estimate the damping 
parameters of vibrations. 

A model of overdamped quantum bath is suggested 
in this paper. It is represented by the spectral den- 
sity, which is directly obtained from the quantum- 
mechanical correlation function of bath coordinates and 
not assuming the classical correlation function as in the 
overdamped semi-classical bath model. The spectral den- 
sity of the overdamped quantum bath decays as uj~ 3 at 
large frequencies, which is preferable to uo~ l used in over- 
damped semi-classical model. In two-dimensional photon 
echo signals, it results in increased homogeneous broad- 
ening of peaks. Dynamics of quantum beats is not af- 
fected by this spectral density; however, the ratio of 
rephasing and non-rephasing signals changes, as well as 
peak intensities. 

Damping of vibrations causes the decay of coher- 
ences; its influence on the peak shape and correla- 



tion/anticorrelation of the diagonal peak width and in- 
tensity is reported to be insignificant compared to the 
description of undamped vibrations. Further increase of 
the damping strength results in instantaneous disappear- 
ance of vibrational coherences and this limit corresponds 
to purely overdamped vibrational motion. 



ACKNOWLEDGMENTS 

This research was partially funded by the European 
Social Fund under the Global grant measure. 



REFERENCES 

1 Y. S. Kim and R. M. Hochstrasser, J. Phys. Chem. B 113, 8231 
(2009). 



10 



2 P. F. Tian, D. Keusters, Y. Suzaki, and W. S. Warren, Science 

300, 1553 (2003). 
3 G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, 

Y. C. Cheng, R. E. Blankenship, and G. R. Fleming, Nature 446, 

782 (2007). 

4 T. Mancal, A. Nemeth, F. Milota, V. Lukes, H. F. Kauffmann, 
and J. Sperling, J. Chem. Phys. 132, 184515 (2010). 

5 D. Abramavicius, V. Butkus, J. Bujokas, and L. Valkunas, 
Chem. Phys. 372, 22 (2010). 

6 H. Okamoto and K. Yoshihara, Chem. Phys. Lett. 177, 568 
(1991). 

7 Y.-C. Cheng and G. R. Flemig, J. Phys. Chem. A 112, 4254 
(2008). 

8 F. Perakis, S. Widmer, and P. Hamm, J. Chem. Phys 134, 204505 
(2011). 

9 E. Collini and G. D. Scholes, Science 323, 5912 (2009). 
10 F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. Plenio, 

J. Chem. Phys. 131, 105106 (2009). 
11 G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. Caram, 

E. Harel, J. Wen, R. Blankenship, and G. Engel, Proc. Natl. 

Acad. Sci. USA 107, 12766 (2010). 
12 E. L. Read, G. S. Engel, T. R. Calhoun, T. Mancal, T. K. Ahn, 

R. E. Blankenship, and G. R. Fleming, Proc. Nat. Acad. Sci. 

USA 104, 14203 (2007). 
13 E. Collini, C. Y. Wong, K. E. Wilk, P. M. G. Curmi, P. Brumer, 

and G. D. Scholes, Nature 463, 644 (2010). 
14 M. Bixon and J. Jortner, J. Chem. Phys. 107, 1470 (1997). 
15 M. H. Vos, J. C. Lambry, S. J. Robles, D. C. Youvan, J. Breton, 

and J. L. Martin, 88, 8885 (1991). 
16 A. Nemeth, F. Milota, T. Mancal, V. Lukes, J. Hauer, H. F. 

Kauffmann, and J. Sperling, J. Chem. Phys. 132, 184514 (2010). 
17 N. Christensson, F. Milota, J. Hauer, J. Sperling, O. Bixner, 

A. Nemeth, and H. F. Kauffmann, J. Phys. Chem. B 115, 5383 

(2011). 

18 D. B. Turner, K. E. Wilk, P. M. G. Curmi, and G. D. Scholes, J. 

Phys. Chem. Lett. 2, 1904 (2011). 
19 V. Butkus, D. Zigmantas, L. Valkunas, and D. Abramavicius, 

ArXiv e-prints (2012), http://arxiv.org/abs/1201.2753. 



20 F. Caycedo-Soler, A. W. Chin, J. Almeida, S. F. Huelga, and 
M. B. Plenio, J. Chem. Phys 136, 155102 (2012). 

21 T. Mancal, N. Christensson, V. Lukes, F. Milota, O. Bixner, 
H. F. Kauffmann, and J. Hauer, J. Phys. Chem. Lett. 3, 1497 
(2012). 

22 N. Christensson, H. F. Kauffmann, T. Pullerits, and T. Mancal, 

J. Phys. Chem. B 0, null (0). 
23 J. M. Womick and A. M. Moran, J. Phys. Chem. B 115, 1347 

(2011). 

24 T. Holstein, Ann. Phys. 8, 325 (1959). 

25 Y.-C. Cheng and R. J. Silbey, J. Chem. Phys. 128, 114713 (2008). 
26 J. Zhu, S. Kais, P. Rebentrost, and A. Aspuru-Guzik, J. Phys. 
Chem. B 115, 1531 (2011). 

27 J. Hu, M. Luo, F. Jiang, R.-X. Xu, and Y. Yan, J. Chem. Phys 
134, 244106 (2011). 

28 J. Prior, A. W. Chin, S. F. Huelga, and M. B. Plenio, Phys. Rev. 
Lett. 105, 050404 (2010). 

29 H. van Amerogen, L. Valkunas, and R. van Grondelle, 
Photosynthetic Excitons (World Scientific, Singapore, 2000). 

30 V. May and O. Kiihn, Charge and Energy Transfer Dynami- 
cs in Molecular Systems (Wiley- VCH, Weinheim, 2011). 

31 S. Mukamel, Principles of Nonlinear Optical Spectrscopy (Ox- 
ford University Press, New York, 1995). 

32 H. -P. Breuer and F. Petruccione, 

The theory of open quantum systems (Oxford University 
Press, New York, 2002). 

33 D. Abramavicius, B. Palmieri, D. V. Voronine, F. Sanda, and 
S. Mukamel, Chem. Rev. 109, 2350 (2009). 

34 V. Butkus, D. Abramavicius, A. Gelzinis, and L. Valkunas, 
Lith. J. Phys. 50, 267 (2010). 

35 D. Abramavicius, L. Valkunas, and S. Mukamel, Europhys. Lett. 
80, 17005 (2007). 

36 V. Butkus, A. Gelzinis, and L. Valkunas, J. Phys. Chem. A 115, 
3876 (2011). 

37 D. Egorova, Chem. Phys. 347, 166 (2008). 

38 A. Gelzinis, D. Abramavicius, and L. Valkunas, Phys. Rev. B 84, 
245430 (2011). 

39 A. V. Pisliakov, T. Mancal, and G. R. Fleming, J. Chem. Phys. 
124, 234505 (2006). 



